Performance Analysis of Queueing Networks via Robust Optimization 

Dimitris Bertsimas * David Gamarnik t Alexander Anatoliy Rikun * 

September 22, 2010 

Abstract 

Performance analysis of queueing networks is one of the most challenging areas of queueing theory. 
Barring very specialized models such as product-form type queueing networks, there exist very few 
results which provide provable non-asymptotic upper and lower bounds on key performance measures. 

In this paper we propose a new performance analysis method, which is based on the robust 
optimization. The basic premise of our approach is as follows: rather than assuming that the 
stochastic primitives of a queueing model satisfy certain probability laws, such as, for example, 
i.i.d. interarrival and service times distributions, we assume that the underlying primitives are 
deterministic and satisfy the implications of such probability laws. These implications take the form 
of simple linear constraints, namely, those motivated by the Law of the Iterated Logarithm (LIL). 
Using this approach we are able to obtain performance bounds on some key performance measures. 
Furthermore, these performance bounds imply similar bounds in the underlying stochastic queueing 
models. 
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Wc demonstrate our approach on two types of queueing networks: a) Tandem Single Class queue- 
ing network and b) Multiclass Single Server queueing network. In both cases, using the proposed 
robust optimization approach, we are able to obtain explicit upper bounds on some steady-state 
performance measures. For example, for the case of TSC system we obtain a bound of the form 
C(l — p)^ 1 lnln((l — /o) -1 ) on the expected steady-state sojourn time, where C is an explicit constant 
and p is the bottleneck traffic intensity. This qualitatively agrees with the correct heavy traffic scaling 
of this performance measure up to the lnln((l — p) ) correction factor. 

1 Introduction 

Performance analysis of queueing networks is one of the most challenging areas of queueing theory. The 
difficulty stems from the presence of network feedback, which introduces a complicated multidimensional 
structure into the stochastic processes underlying the key performance measures. Short of specialized 
cases, such as product form networks, which typically rely on Poisson arrival/exponential service time 
distributional assumptions, the problem is largely unresolved. Specifically, given the topological descrip- 
tion of a queueing network and given the description of the underlying stochastic primitives such as 
interarrival and service times distributions, we do not have good tools for computing exactly or obtaining 
upper and lower bounds on key performance measures, such as, for example average queue lengths and 
waiting times. Some of results which provide non-asymptotic bounds on performance measures can be 
found in jBPT94j . [KK94] . |KM04j . |JOK97j . [BGT96] . [BNM99] . all of which require Markovian (Poisson 
arrival/exponential service time) distributional assumptions. Moreover, some of these bounds become 
quite weak as traffic intensity (of some of the network components) approach unity. For example, a 
bound of the form 0((1 — p*)~ 2 ) is obtained in [BGTOlj . where p* is the bottleneck (real or virtual, see 
the reference) traffic intensity. The other references can lead to infinite upper bounds even in the cases 
where stationary distribution exists. The approaches in these papers also do not extend to the case of 
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non-Markovian systems. As a consequence, most of the known performance analysis results are of an 
asymptotic nature, which apply to queueing networks in various limiting regimes, such as the heavy 
traffic regime [Har90] . [Whi02] . [CYOlj , large deviations methods [GOW04],[SW95j, approximations by 
phase-type distributions [Klc75],[LR87j. 

In this paper, we partially fill this gap by developing a new performance analysis approach based on 
robust optimization methods. The theory of robust optimizaiton emerged recently as a very successful 
and constructive approach for the analysis of certain stochastic modeling problems Soy73j,[BTN98j, 
[BTN99J, |BS03| . |BS04j . The main premise of our approach in the queueing context is that, rather 
than assuming probabilistic laws for the underlying stochastic primitives, such as, for example, i.i.d. 
interarrival and service times, we consider a deterministic queueing model and we will assume only the 
implications of these laws. Specifically we consider implications of the Law of the Iterated Logarithm 
(LIL). The objective is to find laws which on the one hand hold in the underlying stochastic queueing 
model and, on the other hand, lead to linear constraints in the formulation of the robust optimization 
problem, and LIL accomplishes this. We illustrate our approach using two queueing models, namely the 
Tandem Single Class (TSC) queueing system operating under the First-In-First-Out (FIFO) scheduling 
policy, and the Multiclass Single Server (MCSS) queueing system operating under an arbitrary work- 
conserving policy. Motivated by the LIL, we consider constraints of the form ^i<i<fc U% < \~ l k + 
rVHnln/c, for all k > 1. Here (U k ,k > 1) is any of the stochastic primitives of the underlying 
queueing system, such as, for example, the sequence of interarrival times and A stands for the rate of 
this stochastic primitive. Using these bounds, we derive explicit bounds on some performance measures 
such as sojourn time in the TSC system, namely, the time it takes for a job to be processed by all the 
servers, and the virtual workload (virtual waiting time) in the MCSS system, namely, the time required 
to clear the current backlog in the absence of future arrivals. In both models we derive upper bounds 
on the aforementioned performance measures for the corresponding deterministic counterpart models 
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and prove that similar bounds also hold for the same performance measures in the underlying stochastic 
models. In both cases the bounds are of the order 0(j^ In In 377-), where p is the (bottleneck for the 
case of TSC model) traffic intensity. This matches the correct O(j^) order short of lnln((l — p)^ 1 ) 
error. While the technical derivation of these bounds is involved, the conceptual approach is very 
simple. An interesting distinction of our approach from other robust optimization type results is that 
our results are explicit, as opposed to numeric results one typically obtains from the formulating and 
solving a robust optimization model. These explicit bounds however, come at a price of not caring much 
for the constants corresponding to the leading coefficient. In order to keep things simple we sometimes 
use very crude estimates for such constants. 

Our approach bears similarity with some earlier works in the queueing literature. Specifically, the 
pioneering work of Cruz |Cru91a| . |Cru9lb] used a similar non-probabilistic approach to performance 
analysis by deriving bounds based on placing deterministic constraints on the flow of traffic called 
"burstiness constraints". The method could be applied to fairly general network topologies and led 
to more research in the area. In [GP93j,|GP94], tighter performance bounds were obtained assuming 
a "Leaky Bucket" rate admission control from |Tur86| and particular service disciplines. In addition, 
there is some similarity between the philosophy of our approach and the adversarial queueing network 
approach [B KR+01] . [AAF+96] . [Gam03] . [GamOO] . [Goe99j . which emerged in the last decade in the com- 



puter science literature and also replaces the stochastic assumptions with adversarial deterministic ones. 
The deterministic constraints used in the aforementioned works are of the form of A(t) < Xt + B where 
A(t) is the number of external arrivals into the queueing system up to time t and A represents the arrival 
rate. As it turns out, these types of assumptions are too restrictive from the probabilistic point of view 
and do not lead to bounds on the underlying stochastic network: observe that every renewal process 
A(t) arising from an i.i.d. sequence with positive variance violates this assumption almost surely for 
every B for large enough t. As we demonstrate in this paper, the constraints motivated by the LIL, 



namely A(t) < Xt + By/t In hit, can indeed be served to obtain performance bounds, which can be 
translated into the underlying stochastic network. In fact, the key contribution of our approach is that 
the deterministic constraints we place on the service and arrival processes are rich enough to lead to 
stochastic results. The results based on "Leaky Buckets", bounded burstiness and adversarial queueing 
theory address very general queueing networks. It would be an interesting research project to extend 
our results based on robust optimization to these general network structures. 

The rest of the paper is structured as follows. In the following section we describe two queueing 
models under the consideration, namely the tandem single class queueing network and the single server 
multiclass queueing network, as well as their robust optimization counterpart models. Our main results, 
namely the performance bounds in robust optimization type queueing systems and their implications 
for stochastic queueing systems are stated in Section [3J The proofs of our main results are in Sections H] 
and [5j Some concluding thoughts and directions for further research are outlined in Section [6J Several 
technical results necessary for proofs of main theorems are delayed untill the Appendix section. 

We close this section with some notational conventions. In stands for the logarithm with natural 

base. The notation for a non-negative vector x 6 M. d means applying the square root operator 
1 i 

coordinate- wise: (x)a = (a;? , 1 < i < d). A T denotes a transposition operator applied to the matrix A. 

2 Model description 

We now describe the two queueing models analyzed in this paper, both very well studied models in the 
literature. We begin by describing these models in the stochastic setting, and then we describe their 
deterministic robust optimization counterparts. 
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2.1 A tandem single class (TSC) queueing network. Stochastic model 

The model is a tandem of single servers S\,...,Sj processing a single stream of jobs arriving from 
outside and requiring services at Si,... ,Sj in this order. The jobs arrive from outside according to an 
i.i.d. renewal process. Let Ui, U2, U3, . . . denote i.i.d. interarrival times with a common distribution 
function F a (t) = f(U\ < t), where U\ is the time at which the first job arrives. The external arrival 
rate is defined to be A = 1/E[{7i] and the variance of U\ is denoted by a\. 

The jobs arriving externally join the buffer corresponding to server S\ where they are served using 
First-In-First-Out (FIFO) scheduling policy. We assume that all buffers are of infinite capacity. After 
service completion, jobs are routed to the buffer of server S2, where they are also served using FIFO 
scheduling policy, then they are routed to servers S3, S4, etc. After service completion in server Sj the 
jobs depart from the network. Let V£ denote the service time requirement for job k in server j. We 
assume that the sequence (V£,k > 1) is i.i.d. for each j, and is independent from all other random 
variables in the network. The distribution of the service time in server j is F 8 j(t) = F(V± < t),t > 0. 
The service rate in server Sj is defined to be pj = 1/E[V/], and we denote by p m i n = m.mi<j<j p,j 
the rate of the slowest server, cr^ ■ denotes the variance of V( for each j = 1,...,J. The traffic 
intensity in server Sj is defined to be pj = X/fij, and the bottleneck traffic intensity is defined to be 
p* = maxj pj = X/pmin- 

Denote by W 3 k the waiting time experienced by job k in server j not including the service time V£. 
Let Wk = Ylj^Wk + ^fc) ^ e the sojourn time of the job k. Namely, this is time between the arrival 
of job k into buffer 1 and service completion of the same job in buffer J. Denote by Qj(t) the queue 
length in server j (the number of jobs in buffer j) at time t. We assume that initially all queues are 
empty: Qj(0) = 0, 1 < j < J, although most of our results can either be easily adopted to the case 
of non-zero queues at time zero, or apply to the steady-state measures where the initializations of the 
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queues is irrelevant. Let I k denote the idle time of server j in between servicing jobs k — 1 and k for 
k = 2, ...,N. We define l{ = Vj = 1, . . . , J. 

The model just described will be denoted by TSC(St) (Tandem Single Class Stochastic) for short. 
It is known |Sig90 ,[Dai95],[DM95J,|CY0l] that as long as p* < 1, and some additional mild conditions 



hold, such as finiteness of moments, TSC(St) is stable and the stochastic processes underlying the 
performance measures such as queue lengths, workloads, sojourn times are mixing. Namely, these 
processes are positive Harris recurrent [Dai95j,[MT93j, and the transient performance measures converge 
to the (unique) steady-state performance measures both in distributions and in moments. Computing 
these performance measures is a different matter, however. We denote by Woo,Woa the steady state 
versions of the random variables W° k , Wf, . Thus provided that p* < 1 and some additional technical 
assumptions hold, we have 

lim E[W n ] = E[Woo]. (1) 

n— >oo 

We will assume that p* < 1 holds without explicitly stating it. Rather than describing the assumptions 
required to make ([T|) true, we will simply assume when stating our results that (pQ) holds as well. 

2.2 A multiclass single server (MCSS) queueing system. Stochastic model 

We now describe our second queueing model. Consider a single server queueing system which processes 
J classes of jobs. The jobs of class j = 1, 2, . . . , J arrive from outside according to a renewal process 
with i.i.d. interarrival times Ul,k > 1 and distribution function F a j(t) = ¥(U( < t). The arrival 
rate for class j jobs is Xj = 1/K[U(}. It is possible that some classes j do not have an external arrival 
process, in which case U 3 k = oo almost surely and Xj = 0. Let a^j be the variance of U{. The sequences 
(U k ,k > 1) are also assumed to be independent for different j. Let A = (Xj) denote the J-vector of 
arrival rates. We let A max = maxi< 3 < j Xj and A m i n = mm i<j<jAj. We let A(t) = (Aj(t)) denote the 
vector of cumulative number of external arrivals up to time t where Aj(t) = max{/c : ^i<j<fc < t}. 



The jobs corresponding to class j are stored in buffer Bj until served. As in the single class case, 
we assume all buffers are of infinite capacity. The service time for the k-th job arriving to buffer Bj 
is denoted by V£ and the sequence (Vj!,k > 1) is assumed to be i.i.d. with a common distribution 
function F s j(t) = P(V^ < t). Additionally, these sequences are assumed to be independent for all j 
and independent from the interarrival times sequences {U^,k > 1). The average service time for class 
j is rrij = ^[Vj 7 ] and the service rate is fAj = l/EfVj 3 ]. o~ 2 s • denotes the variance of V(. Let m = {rrij) 
denote the J- vector of average service times and let /i = (fj,j) be the J- vector of service rates. Let M 
denote the diagonal matrix with j-th entry equal to fj,j and let ^ max = maxi<j<j fij. 

We assume that the jobs in buffer Bj are served using FIFO rule, but prioritizing jobs between 
different buffers Bj is done using some scheduling policy 6. The only assumption we make about 8 is 
that it is a work- conserving policy. Namely, the server is working full time as long as there is at least one 
job in one of the buffers Bj, 1 < j < J. The only performance measure we will consider is the workload 
(defined below) for which it is well known that the details of the scheduling policy are unimportant for 
us, as long as the policy is work-conserving. 

The routing of jobs after service completions is determined using a routing matrix P, which is an J 
by J 0, 1 matrix P = (Pij, 1 < i,j < J). It is assumed that Y2j Pi,j — 1 f° r each i. (Namely, the sum is 
either 1 or 0). Upon service completion in buffer Bi, the job of class i is routed to buffer j if Pij = 1. 
Otherwise, if ^ . Pij = 0, the jobs of class i leave the network. It is assumed that P n = for some 
positive integer n. It is easy to see that this condition is equivalent to saying that all jobs eventually 
leave the network. 

It is known [CYOlj that the traffic equation \ = Aj + X^i<j<j nas a unique solution A = (Xj) 

given simply as A = [I — P T ] _1 A, where I is the J by J identity matrix. Let A max = maxj(Aj) (observe 
that Xj < Xj for every j and hence A max > A max ). Let A(t) = (Aj(i)) denote the vector of number of 
arrivals by time t that will eventually route to server j: Aj(t) = + (P T ) 1 + (P T ) 2 + • • -)A(t) = 
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ej[l — P T ]~ 1 A(t) and ej denotes the j — th unit vector. 

The traffic intensity vector is defined to be p = M -1 A = M _1 [J — P T ]~ 1 X. The traffic intensity of 
the entire server is p = e T p, where e is the J vector of ones. Let Qj(t) denote the queue length in buffer 
j at time t, let Q(t) = (Qj(t)). We assume that Q(0) = 1. As for the case of TSC model, our results can 
be extended to the case Q(0) > 0, but for the results regarding steady-state behavior, the initialization 
of queues is irrelevant. Denote by Wl the waiting time of the fc-th job arriving into buffer j. We let 
Wf denote the workload at time t. Namely, Wt is the time required to process all the jobs present in 
the system at time t, in the absence of the future arrivals. Note that Wt is also the virtual waiting time 
at time t when the scheduling policy is FIFO. Observe that if to marks the beginning of a busy period 
and t\ belongs to the same busy period (namely, the server was working continuously during the time 
interval [io,ii]), then almost surely 

w tl = £ vi + ...+ Yl yj-ih-tv). (2) 

i=Ai(to) i=Aj(t ) 
The model described above is denoted by MCSS(St) (Multiclass Single Server Stochastic) for short. It is 

known [Dai95j that if p < 1, and some additional technical assumption on inter arrival and service time 

distributions hold then MCSS(St) is stable and enters the steady state in the same sense as described 

for the tandem queueing network. While in this case the steady-state distribution of many performance 

measures usually depends on the details of work-conserving policy used, the steady-state distribution of 

the workload does not depend on the policy, as we have discussed above. Let denote the workload 

in steady state, and let and 1^ denote the steady-state duration of the busy and idle periods, 

respectively. Additionally, denote by Jo, B\, I\, E>2, 12, ■ ■ ■ the alternating sequence of the lengths of the 

busy and idle periods of the MCSS(St) system, assuming that time zero initiates a busy period. Under 
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the same technical assumptions as above the following ergodic properties hold almost surely: 

lha I^l =E[Wool (3) 

t— >oo t 

lim El ^ nBi =E[B 00 ], (4) 

lim El ^ / * =E[/ 0O ], (5) 

lim El <*<"^ = E[£gJ. (6) 

n— >oo tt, 

We denote by n(t) the number of busy periods that have been initiated up to time t. Mathematically, 
we define n(t) to satisfy T,i<i< n (t)-i( B i + J i) < * < Ei<i<„(t) O 8 * + T i)- When 1 6 Ei<i<n(t)-l( B i + 
^)>Ei<i<n(t)-l(^» + ^) + -^n(t)]) * f &ns on a busy period and using the definition of n(t), we have 
< B n{t) . When t E [£ 

l<i<n(t)-l(-^» + + El<i<ti(t)(^i + * ^ ans on idle P er i°d / n (t) 

and hence W(t) = 0. We let r, denote the beginning of the i-th busy period. This implies 

£ w(s)ds _ e£? £f** +fl «'*> < Ex^t) ^ 2 



* * El< 1 ;<n(()-l( B i + 7 ») 

If ©,©,© and © hold, then we also obtain 

HWoolS nB ao ]+E[I 00 ] S E[B oo y {V 
This bound will turn useful when we apply our results for robust optimization models to the underlying 
stochastic model. As for the TSC case, we assume from now on p < 1. Rather than listing the 
assumptions leading to ergodic properties ((3j) , (S]) , dS]) and ([6]) we assume when stating our results, that 
the stochastic process Wt enters the steady-state as t — > oo and that the properties j3]),(J2|), ([5]) and ([6]) 
holds almost surely. 

2.3 Robust optimization type queueing systems 

We now describe deterministic robust optimization type counterparts of the two stochastic queueing 
models described in the previous subsections. 
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We begin with TSC model and describe the corresponding model which we denote by TSC(RO) 
(Tandem Single Class Robust Optimization). The description of the network topology is the same as 
for TSC(St). However, it is not assumed that U^, Vj! and, as a result Q(t), W%, are random variables. 
Rather we assume that these quantities are arbitrary subject to certain linear constraints detailed below. 
Additionally, we assume that the system starts empty Q(0) = and only n jobs go through the system. 

Specifically, consider a sequence of non- negative deterministic interarrival and service times (Uk, 1 < 
k < n), (V/}, 1 < k < n), 1 < j < J. Let 



(x)={ 



s/x\n lnx, x > e e ; 

(8) 

1, x < e e . 



We assume that there exist A, T a and fJLj,T s j > 0, 1 < j < J such that 

| U k~ A _1 (n- k)\ <T a (j)(n- k), k = 0, 1, . . . , n - 1, (9) 

k+l<i<n 

| Yl Vi-pq\n-k)\ <T s>j( f>(n-k), k = 0, 1, . . . , n - 1, j = 1, 2, . . . , J. (10) 

k+l<i<n 

It is because we need to consider tail summation Ylk+Ki<n we assume that only n jobs going through 
the system, though we will be able to apply our results in the stochastic setting where infinite number of 
jobs pass through the system. Let T = max(r a , T s j). Borrowing from the robust optimization literature 
terminology ([B S04| ). the parameters r a ,r sj ,r are called budgets of uncertainty. Note, that the values 
Uk,Vj!,k > 1 uniquely define the corresponding performance measures Qj(t),W£,Wk,k = l,...,n. 
There is no notion of steady state quantities Qj(oo),W OQ for the model TSC(RO). The motivation for 
constraints (|9|) and (jlpp comes from the Law of the Iterated Logarithm, and we discuss the connection 
in a separate subsection. 

We denote the robust optimization counterpart of the MCSS(St) model by MCSS(RO). In this case 
it turns out to be convenient to consider infinite sequence of jobs. Thus consider infinite sequences 
of deterministic non-negative values (Ul,k > 1), (Vi,k > 1),1 < j < J. It is assumed that values 
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Xj, fj,j, T a j, r s j > 0, 1 < j < J exist such that 

\J2 U k~ X I %k \ ^ r «J<t>( k )> k = l,2,... , j = 1, 2, ... , J, (11) 

l<i<fc 

\Yj V i~ h lk \ ^ r «^( fc )> k = 1,2, . . . , j = 1,2, . . . , J. (12) 

l<i<fe 

For convenience we assume that at time zero the system begins with exactly one job in every class 
j = 1, . . . , J: Qj(0) = 1. Then the first after time zero external arrival into buffer j occurs at time U(. 
As before, we let T = max(r aj -, T s j). 

For technical reasons, we also assume that T in TSC(RO), MCSS(RO) constraints satisfies 

Ar > e 2e and minAjT>e 2e , respectively. (13) 

j 

2.4 The Law of the Iterated Logarithm 

One of the cornerstones of the probability theory is the Law of the Iterated Logarithm (LIL) |Chu01| . 
which states that given a i.i.d. sequence of random variables Xi, . . . , X n , . . . with zero mean and finite 
variance a, the following holds almost surely, 

,. J2l<k<n^-k , Yll<k<n^-k 

iimsup — ~ = = = 1, hmmf — ~ = = = —1. 
n->oo crv2nlnlnn n->oo aV^nlnlnn 

The LIL extends immediately to non-zero mean i.i.d. sequences by subtracting raELYi] from J2i<k<n -^k- 
Furthermore, LIL implies (in the case of zero-mean variables) that 

-p A \J2l<k<n X k\ 

Tlil = sup - < co, (14) 

n>i ay/ 2(j){n) 

where (j) is defined in ([8]). Note that Tlil is a random variable. Thus when we consider stochastic 
queueing models such as TSC(St) or MCSS(St), the constraints ©, CEROID, 02} hold with probabil- 
ity one, with T = ^/^.T^i^a, where Tlil is defined in (|14p for the corresponding random sequence. 
Specifically, let T a = T^lil = Tlil and a = a a , when Xk = U n -k — A _1 ,0 < k < n — 1 and Uk 
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is the sequence of interarrival times in the TSC(St) model. Similarly define T s j = T s jLlL when 
Xk = Vn-k ~ ^J 1 '® — ^ — 71 ~ ^ — 3 — J' Observe, that for T a ,T s> j thus defined, the constraints 
([9]). (|10p hold for an infinite sequences of jobs (that is jobs which would have indices — 1, — 2, . . .), even 
though we need it only for the first n jobs. For the MCSS(St) model define T a j = r ajj LiL, r sj - = r sjj LiL 
corresponding to the sequences U J k — XJ , V£ — XJ , k > 1, respectively. We obtain 

Proposition 1. Constraints ^) , [W \) . [U\) . [T3 \) hold with probability one for T a = V2T ' a ,LiLO~a, r sj - = 
V^ s ,j,LiLO- s j, T a j = \/2T a j iLI LO- a ,j, andT s j = y/2T s ,j,LiL<Ts,j, respectively, where T.^. yLIL is defined in 
for the corresponding sequence. 

As a conclusion, for every property derivable on the basis of these constraints in our deterministic 
robust optimization queueing network models, such as, for example, bounds on the sojourn time of the 
n-th job in TSC, the same property applies with probability one for the underlying stochastic network. 
This observation underlies the main idea of the paper. 

3 Main results 

In this section we state our main results on the performance bounds for robust optimization type 
queueing networks TSC(RO) and MCSS(RO), and the implications of our results for their stochastic 
counterparts TSC(St) and MCSS(St). We begin with TSC(RO) with the goal of obtaining a bound on 
the sojourn time. 

Theorem 1. The sojourn time of the n-thjob in the TSC(RO) queueing system with constraints IPj) A1U\) 
satisfies 

W n < 7 -f^lnln-^ + JX-\ (15) 
1 — p* 1 — p* 

Observe that the bound on the sojourn time is explicit. It is expressed directly in terms of the 
primitives of the queueing system such as arrival and service rates. Observe also that the upper bound 
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is independent from n. One can think of this bound as a "steady-state" bound on the sojourn time in the 
robust optimization model of the TSC system. Additionally, the constant T 2 is related to the "variances" 
of interarrival and service times viz a vi the LIL (|14p . It is known that in the stochastic GI/GI/1 queueing 
system the expected waiting time in steady state is approximately (<7 2 +<7 2 )/(2A(l— p)), when the system 
is in heavy traffic, namely p — > 1. Namely, the expected waiting time depends linearly on the variances 
of interarrival and service time. Our bound ([To]) is thus consistent with this type of dependence. On 
the other hand, unfortunately, our bound depends quadratically on the number of servers J, whereas 
the correct dependence is known to be linear, at least in some special cases [Rei84j,[GZ06j. 

The bound above does not have a correct 0((1 — scaling, which is known to be correct from 

the heavy-traffic theory perspective [ Rei84 j . [GZ06] . However, the correction factor is a very slowly 
growing function In In. The upshot is that we can use this bound to obtain a bound on W n and Woo in 
the underlying stochastic system. This is what we do next. 



Corollary 1. For every n > 1 the sojourn time of the n-thjob in the TSC(St) queueing network satisfies 

+ J\- 1 . (16) 



E[W n \ < E 



7J 2 r 2 A, , jxf 

In In ■ 



1 — p* 1 — p* 

where F = maXj{^\/2a a F a lil, V^&s jF s ,j lil, e 2e \~ l ). If in addition the assumption |7]) holds th 



en 



IE WJ < E 



7J 2 r 2 A, , JXF 
In In 



+ J\- L . (17) 



1- p* 1 - p* 

Proof. We first assume Theorem [T] is established. Note, in the context of the stochastic system, both 
W n and T in Theorem [T] are random variables. We take F = maxj(\/2cr a r aj LiL, V^^jT^lil, e 2e A _1 ) 
to satisfy (fT5|) . where r V) LiL is defined in (fT4"|) for the corresponding sequence. Applying Proposition Q] 
we have that (fl~5|) holds with probability one for the underlying stochastic network. The bound (fT6|) 
now follows from taking expectations of both sides of (I15p . The bound (|17p follows from applying ([TJ 
to (US]). □ 
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We now turn our attention to the MCSS queueing model. Our approach for deriving a bound on the 
workload is based on first obtaining an upper bound on the duration of the busy period. Thus, we first 
give a bound on the duration of the busy period and then turn to the workload. Recall our assumption 
Q(0) = 1, though our results can readily be extended to the general case of Q(0) > 0. Thus, time t = 
marks the beginning of a busy period. 

Theorem 2. Given a MCSS(RO) queueing system with constraints (G2P, (dp? let B be the duration of 
the busy period initiated at time 0. Then 

K^^yhh""^, us) 

(i - py i - p 

and sup W(t) < 2 ( 4J + 3 ) 2A ma x r 4 lnln (4J + 3)AL x r 2 + p + ^2 T 3 (19) 

0<t<B 1 - P 1 - P 

While the bound (I19p corresponds to the maximum workload during a given busy period, the actual 
value of the bound does not depend on the busy period length explicitly. As it will become apparent 
from the proof, we use the same technique for obtaining a bound simultaneously on the duration of the 
busy period and maximum workload during the busy period. Let us now discuss the implications of 
these bounds for the underlying stochastic model MCSS (St). 

Corollary 2. Given a MCSS(St) model, suppose the relations (0), (E|) and {6}) hold. Then 

E[Boo] < E rM4J + 3)*ALj' ]n ln 2(47 + 3)^1 
L (1-p) 2 1-P J 



25(4 J + 3) 4 ALx/Uma x rV lnln 2(4 J + 3)^^x2 



(21) 



E ^ E L (!-,)« V 1-p 

where T = m&Xj(V2o- a ,jT atjjLIL , \/2a s jT ' s ,j,LiL, e 2e A~[ n ). 

Unfortunately, in this case the scaling of our bounds as p — >• 1 deviates significantly from the 
correct behavior. Prom the heavy traffic theory [DK95J, the correct behavior for the steady-state 
workload should be 0((1 — p) -1 )- As for the steady-state busy period, the theory of M/G/l queueing 
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system |Kle75| suggests the behavior 0((1 — /o)~ 5 ) as opposed to 0((1 — /5)~ 2 lnln(l — p)~ l ) which 
we obtain. On the positive side, however, we managed to obtain explicit bounds on the performance 
measures which are expressed directly in terms of the stochastic primitives of the model, which we do 
not believe was possible using prior methods. We leave it as an interesting open problem to derive 
the performance bounds based on the robust optimization technique, which lead to the correct scaling 
behavior as p — > 1. 

While the proofs of our main results are technically involved, conceptually they are not complicated. 
Before we turn to formal proofs, in order to help the reader, we outline below informally some of the 
key proof steps for our results. 

For the TSC queueing network we first replace the constraints (f9|). (fTUj) with more general constraints, 
see (j22|) and (|23p below. Our results for the TSC network rely mostly on the Lindley's type recursion 
which in a single server queueing system recursively represents in the waiting time of the n-th job in 
terms of the interarrival and service times of the first n jobs. It is classical result of the queueing 
theory that this waiting time can be thought of as maximum of a random walk, with steps equalling in 
distribution to the difference between the interarrival and service times. We derive a similar relation 
in the form of a bound on the sojourn time of the n-th job in the TSC network. This bound is given 
in Theorem [3J Then we view this bound as an optimization problem and obtain a bound on the 
objective value by proving the concavity of the objective function and substituting explicit bounds from 
constraints (f9lh (fT0|) . 

Our proofs for the MCSS queueing system rely on the relation (J2J) . Namely, we take advantage of the 
fact that the workload is depleted with the unit rate during the busy period. Then we take advantage of 
the constraints f)ll[) . fjl2|) to show that in the MCSS(RO) system the workload at time t during the busy 
period can be upper bounded by an expression of the form —at + bV t In In t + c with strictly positive a, b. 
It is then not hard to obtain an explicit estimated to such that this expression is negative for t > to- 
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Since this expression is an upper bound on a non- negative quantity (workload), then the duration of 
the busy period cannot be larger than to. This leads to an upper bound on the duration of the busy 
period in the MCSS(RO) system. In order to obtain a bound on the workload, we again take advantage 
of ([2]) and further obtain explicit upper bounds on the terms involving the sums of service times. We 
show that the workload at time t is at most — at + bV t In In t + c. We then obtain an upper bound on 
the workload during the busy period by obtaining explicit bounds on max t > — at + &V t In In t + c. 

Our derivation of the bounds for the stochastic model MCSS(St) relies on the ergodic representation 
(|3|). We consider a modified system in which each busy period is initiated with simultaneous arrival 
of one job into every buffer j. This leads to a alternating renewal process with alternating i.i.d. busy 
and idle periods. We then obtain a bound on the steady-state workload in terms of the second moment 
of the busy period in the modified queueing system, using the renewal theory type arguments. It is 
this necessity to look at the second moment of the busy period which leads to a conservative scaling 

in our bound (|2ip on the steady-state workload. 

4 Tandem single class queueing system analysis: proof of Theorem [I] 

In order to prove Theorem Q] we first generalize constraints ([9|). (|10p and obtain a method for bounding 
W n under more general uncertainty assumptions. 

4.1 General upper bound on the sojourn times 

Given a sequence of non-negative real values Y 3 m - m [k), rjqax(fc) 1 < j < J, 1 < k < n, r mm (/c), r max (/c) 
1 < k < n, we consider the set of all sequences of service times and interarrival times (V/), (Ui) 
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j = 1, . . . , J, i = 1, . . . , n satisfying for all k = 1, . . . , n 

n 

r J mi „W<E^ r max(£0, (22) 
i=k 
n 

r mi „(£;)<^^<r max (£;), (23) 

i=k 

V?,Ui>0. 

In the next theorem we obtain a bound on the sojourn time of the n-th job in TSC(RO) system in 
terms of values r^ in (fc), T J max (k), T min (k), T max (A;). 

Theorem 3. Suppose the relations Iff)) and i ffgj) hold. Then 

J-l 

W n < max J2 ( r 'max(%) - r J min (%+i + 1)) + rLx(^) - IWA* + 1) (24) 

n>kj>...>ki>l £ — ' 

We now show how Theorem [3] implies our main result Theorem [TJ 

Proof of Theorem [7J The proof consists of two steps: the first step uses Theorem to bound W n with 

uncertainty sets ([9|). (|10p . The second step involves solving some associated maximization problem. 

We set r min (A;) = \~ l (n+l-k)-r a (t>(n + 1 - £;),r max (£;) = A" V+l-A;)+r a (£(n + 1 - fc), r J min (fc) = 

/ij 1 (n + 1 — A;) — T s j(p(n + 1 — k), r J m£K (fc) = /i~ 1 (n + 1 — fc) + T s j(f)(n + 1 — k), where is defined by 

([8]). From Theorem [3] we obtain: 

J-i 

W„ < max V (fi^On + 1 - fc<) + T s Mn + 1 - kj)) 

n>kj>. .. >fci>l ^ V J J / 

J-l 

- 2 + 1 " ~ ^ ~ r ^( n + 1 " " ^ 

+ (^(n + 1 - fej) + r SjJ 0(n + 1 - kj)) - {\-\n + 1 - h - 1) - r o 0(n + 1 - h - 1)) 
Since n > kj + i > kj Vj, we can replace nj 1 by ^ m ; n = max(/i7 , /«2 1 • • • ^j 1 ) < A -1 and preserve 
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inequality. Similarly, we can replace r S) i, r S) 2, . . . , T 8) j, T a by T. We obtain: 

J-l 

W n < max Yl knin(^'+i + 1 " + ^{n + 1 - kj) + </>{n - k j+ i)) 

n>kj>...>ki>l — L 
J-l 

+ (/Vi> + 1 ~ fcj) + I>(n + 1 - fcj)) - (A- X (n - fca) - I>(n - h)) 



< max M~- n (n - h) + 2JT<j>(n + l-k x ) 

n>k\>l 



+ <//i m ; n — A 1 (n — fci) where we used k\ < k% < . . . < kj to combine T terms 



max (n + 1 - - A ) + 2JT(f>(n + 1 — ki) + (J — l)p mn + A 

n>fei>l 



< max (n + 1 - AiiX/C^ - A" 1 ) + 2jr0(ra + 1 - A*) + JA" 1 since A" 1 > jC? 

n>fci>l 



We let x = n + 1 — fei. Since 1 < k\ < n we have that 1 < x < n and obtain: 



W n < max — A ) + 2JT(f ) (x) + JA 

n>x>l 



< max x(//l] n - A" 1 ) + 2JT^(x) + JA -1 (25) 

X>1 

Putting a = A -1 — p~^ in , £> = JT, c = JA -1 , and using the assumption (fT3|) . we have b/a = AJT / (1— p*) > 
e 2e , namely, the condition (|34p is satisfied. Applying Proposition [2] from Appendix we obtain 

7AJ 2 r 2 , , AJr „ , 

W n < In In h JA" 1 . 

1- p* 1 - p* 

This completes the proof of the theorem. □ 
4.2 Proof of Theorem [3] 

Job 1 enters the system first, followed by jobs 2,3, ... ,n. Let U\ be the time between the arrival of 
job i and job i — 1 into server j for j = 2, . . . , n and j = 1, . . . , J. Specifically, U\ = Ui, and we define 
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Ul =V( for j = 2, . . . , J. The following relations are well known in the queueing theory |Kle75| . 

W{ = max(Wti + ^ii -Ui,0) V i = 2, . . . ,n, j = 1, . . . , J, (26) 

/■/ 1 • // ' Vt = 2,...,n,j = 2,...,J, (27) 

Wf = max{ max ^(^'-C^ +1 ),0} V i = 2, . . . ,n, j = 1, . . . , J, (28) 

1 l=k 

liy ; liy // il;', /./) V* = 2,...,n, j = l,...,J. (29) 

We now prove some more detailed results regarding the dynamics of our queueing system. 
Corollary 3. The following relations hold for k = 2, . . . , n — 1: 

E u i = E ctf+tf) = w2-w2+ e ^ + 

i=fc+l i=fc+l i=fc+l 

Proof. The first equality follows from ()27p . To prove the second equality we use (|29p to obtain 

n n 

E^ 1+/ ')= E (^-^-i+^+k 1 -^ 

i=k+l i=k+l 

and the result follows. □ 
Lemma 1. 

k2 k-j kj n n 

= n>t ™* ti>l E 1/ i ' + E l '. 2 + ---+ E tf- + E ^ - E < 3 »> 

i=k\ i=k2 i=kj_i i=kj j=fci + l 

Proof. We prove Lemma Q] by induction. We let W( ,S = + denote the sojourn time of customer 
i in server j. 

Case J = 1: We first define Yli=j+i = f° r an i- Using (|28l) and V^ J > we have for any n = 2, . . . , n: 

n-l 

>fciL- 

i=k 

i 

c 5 

>fci>i 

i=ki j=fci + l 



i=fc 

ra n 

i=k± j=fci + l 

n n 

^ V^ 1 — ^2 and this completes case J = 1. 



max 

n>fci>_ 

i=fci i=fci+l 
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Case J > 1: Note that W n = W n ' + (W n ' + . . . + W n ' ) and denotes the sojourn time of job n in 
J-server system. We suppose that the result holds for a J — 1 tandem system and proceed by induction: 



k2 k 3 

max 

n>k 



rj 
• i 



-- 1 . 1>1 E". 1 + Ek 2 + ...+ E tf-' + Evf- E vi 

i=k\ i=k2 i=kj_i i=kj i=k\+l 

fc 2 k 2 n k 3 kj 

^^(E*?- E E d? + 5: »?+...+ E r + E^ 

j=fci i=fci+l i=k 2 +l i=k 2 i=fej_i i=&j 

fc 2 A; 2 n fc3 fcj n 

-sUE". 1 - E «?)- E c/ + E". 2 +-+ E tf-' + Ei 

i=fel i=fci+l i=fc 2 +l i=k2 i=kj-i i=kj 

n kg kj n 

max Wl' s - V U} + V Vf + . . . + V V/' 1 + V K J the base case J = 1 is used 

n>kj>...>k 2 >l k2 ^ ^ ^ % ^ 

i=k2+l i=k2 i=kj-i i=kj 

n k-j, kj n 

,s iai K s ~ E + E ^ + • • • + E vf- 1 + E k j 



max 

n>fcj>...>fc 2 >l 



n>k 



i=k 2 +l i=k 2 i=kj_i i=kj 

n k 3 kjn 



.1 



- ,> l S,>,(*-' s - E «f) + E »?+■■■+ E tf-' + E* 

i=k 2 +l i=k2 i=fej_i i=kj 

we used Corollary [3] and = + VjJ 

&3 fej n n 

= max Vi; 2 • ... • V F/^+VF/- V c/ 2 

i=fc 2 i—kj— i i=kj i=fc 2 +l 

= W 1 ' + (W n ' + . . . + W n ' ) by inductive assumption on J — 1 server system 
and the proof follows from definition of sojourn time W n . □ 
Proof of Theorem The result follows immediately from Lemma [TJ □ 

5 Multiclass single server analysis: proofs of main results 

5.1 Proof of Theorem [2] 
Lemma 2. For every t satisfying 

t > max(A7V, AT 1 + SX^X^T 2 ), (31) 
the following holds: Aj(t) < tXj + 3X 2 T 2 4>(tXj). 
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Proof. Assume first Aj(t) < e e . Then applying (jlip corresponding to the case Aj(t) < e e , we obtain 
^i( t ) A 7 1 - T a,j < t, namely Aj(t) < Xjt + XjT a j < Xjt + A-,T. Since XjT, 4>(tXj) > 1 from fFjft and ©, 
the desired result is obtained. For the rest of the proof assume Aj(t) > e e . Applying (fTT]) . we obtain 
Aj(t)X~ l - Tajy/Ajit) In In A,- (t) < t. Which gives 



Aj(t) - tXj 
y/Aj(t) \n\n Aj(t) 



<XjT aJ <XjT. (32) 



Define bj by: bj = tXj + 3A 2 r 2 y/tXj In In t Xj . Observe that: 

bj - tXj 3Afr 2 ya7h7ElA7 



V 6 i lnln6 i ({tXj + 3X 2 T 2 y /tX j ]n\ntX j ) \nln(tXj + 3X 2 j r 2 y /tX j \n\ntX j )^ 2 



3X?T 2 y/tXj]n\ntXj 

> j- since iAj > lnlntAj for tXj > e e from 



(tAj + 3A|ryt2 A 2) m l n ( tAj + 3A 2 r2 J t 2\2) 



3A 2 r 2 v / U7lnlniA7 



(tA,)(l + 3A 2 r 2 ) lnln (tXj)(l + 3A 2 r 



3A 2 rViA ? -lnln£A,- 
> ^ — r since tXj > 1 + 3A 2 r 2 from §31 

(tXj)(l + 3A 2 r 2 )lnln(iA i ) 2N " 



3A 2 r 2 /lnln^A 7 • 
> , since 2 lnlntAj > lnln(tAj) for tXj > e e and XjT > 1 

J(4A 2 r 2 )(21nlntA j ) 
> AjT by simplifying above expression. 

Since f ^ is an increasing function for x > e e and from (1321) . we have that > A,- ft) and the result 

Va; In In a; ~~ * ' J — ■> x ' 

is obtained. □ 

We now obtain an upper bound on the cumulative arrival processes Aj(t), 1 < j < J. 
Lemma 3. For every t satisfying [31\) . the following holds 

0(A J (t))<((2 + 6A 2 nax r 2 ))^(A,t) 

Proof. Consider first the case Aj(t) < e e . From (jHJ), we have that (j)(Aj(t)) = 1 and applying (f3Tj) . the 
lemma follows. Now we consider the case Aj(t) > e e . Recall that Aj(t) = ej[I - P T ]- l A(t). Applying 
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Lemma [2] 



Aj(t) < ej [I - P 1 ]- x At + ej [I - P 



7T1-I 



3A?r 2 <^(tAi; 

3Alr 2 ^(tA 2 ^ 



3A2r 2 0(Uj) 



< ej[l - P T ] -1 At + 3A 2 nax r 2 eJ[J - P T ] _1 At applying Ijgip and x > <j>(x) for x > 
= Xjt(l + SA^^r 2 ), applying the definition of \j. 

Applying this bound we also obtain 

lnlnl,(i)< mln^l + 3X 2 max T 2 )) 

< lnln(Ajt) 2 using assumption (f3Tj) 
= In In A j t + In 2 

< 2 In In A^, using Xjt > Xjt > e e from (|3T|) . 
Combining the previous bounds with definition of 4>(x), the lemma follows. 

Lemma 4. For every t satisfying (E2P, we have: m T A(t) — t < (p — l)t + 3A max r 2 (/>(A max i). 
Proof. Applying definition of Aj (t) , we have 

m T A{t) - t = fn T [I - P T ]~ l A{t) - t 

< m T [I - P T ]' 1 (Xt + 3A max r 2 </>(A max t)A) - t from Lemma [2] 

= m j^jt + 3A max r 2 0(A max t) m j^j ~ * applying the definition of A 

3 i 

= {p- l)t + 3A max r 2 </)(A max t) / 9 

and the lemma follows from applying the condition p < 1 to the second term. 

23 



We now obtain an upper bound in the duration of the busy period. Recall the identity ([2]). Since 

the busy period begins at time zero its duration is upper bounded by the first time t such that 

Ai(t) Aj(t) 

J2v' + ...+ £ V/-t<0. (33) 
i=i i=i 

Consider any t satisfying the lower bound (]3ip . We have 
Ai(t) Aj(t) 

i=l i=l 

J J 

< ^ Mj^iC*) + Yl ^MM-t)) - t applying (HU),® 

3=1 3=1 

< m T A(i) - i + T aj ((2 + 6A^ ax r 2 )) 2 0(Ajt) applying Lemma [3] 

3=1 

J 

< t{p - 1) + 3A max r 2 </»(A max t) + J^r(2 + 6A 2 iax r 2 )5(/>(A i t) applying Lemma [1 

3=1 

< t(p - 1) + (4 J + 3)A max rV(A max t), 

where we have used a crude estimate 2+6A 2 nax r 2 < lGA^^r 2 , justified by (fl"3j) . We now apply Lemma[6] 
with x = A max i, a = A~ ax (l — p),b = (4J + 3)A max T 2 /2 and c = 0. The condition (J33J) is implied by 
assumption (fl~3]). and the second condition of Lemma [6] is satisfied since c = 0. We obtain that (|33|) 
holds for all t satisfying (f3Tj) and 

t > 18(4 J + 3) 2 AL x r 4 lnln 3(4 J + 3)A max r 2 

4 A max Amax (1-p) 2 2AmL(l-p) 

5(4 J + 3) 2 Al ax r 4 , , 2(4J + 3)AL x r 2 
{1-pf 1-p 

Observe using (|13p that the right-hand side of the last expression is larger than the right-hand side of 
(f3Tj) . Combining two cases we obtain (fl~8|) . 

We now turn to (fTUj) . First suppose t does not satisfy (fBTjh Denote the right-hand side of (f3"Tj) by 
C. That is t < C. Observe that W(t) < (C - t) + W(C) < C + W(C) as the workload at time C 
corresponds in addition to arrivals during [t, C]. So now we focus on the case when t satisfies (131j) . We 
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use Proposition [2] from Appendix and obtain 



C<t<B 



sup W(t) < 



< 



7(4J + 3) 2 A^r 4 

4AmL(l - p) 

2(4J + 3) 2 AL x r 4 
1-p 



In In 



In In 



(4J + 3)A max r 2 

2AmL(l - p) 
(4J + 3)AL x r 2 



1-p 



From (fl"3j) . we have V > A min . We conclude that 




In In 



(4J + 3)Al ax r 2 
1-p 




This completes the proof of the theorem. 
5.2 Proof of Corollary H 

First we establish bound (|20p . Let t = mark the beginning of a busy period with (random) length 
Boo in steady state. This means that there is an arrival into one of the classes jo at time 0. Consider 
a modified system where the first arrivals into classes j ^ jo, Xj > after time are artificially pushed 
down to exactly time 0. Namely, now at time zero there is an arrival into every class j with Aj > 0. 
The subsequent arrivals into these classes are also pushed earlier by the same amount, thus creating an 
i.i.d. renewal process initiated at time 0. Let B be the busy period initiated in the modified system at 
time 0. It is easy to see that almost surely B > Boo. However, now that we have arrivals in every class 
at time zero, applying Proposition Q] and our result for the robust optimization counterpart queueing 
system, namely applying part (fTBj) of Theorem [21 we obtain the required bound by taking the expected 
values of both sides of (fTHj) . This establishes part ([20|h 

In order to prove (f2T|) . we use a bound ([7]). Using our earlier argument for the proof of (f20|) but 
applying it to the second moment of B we obtain 



On the other hand, we trivially have have E[5oo] > mini<j<j rrij = l/p m ax, since every busy period 
involves at least one service completion. The result then follows. 



E[B£,] < E[f? 2 ] < E 



25(4J + 3) 4 AL x r 8 
(1 - P) 4 



( 



In In 



2(4J + 3)A 2 nax rV 
1-p J 
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6 Conclusion 

Using ideas from the robust optimization theory we have developed a new method for conducting per- 
formance analysis of queueing networks. The essence of our approach is replacing stochastic primitives 
of the underlying queueing system with deterministic quantities which satisfy the implications of some 
probability laws. These implications take the form of linear constraints and for the case of two queueing 
systems, namely Tandem Single Class queueing networks and Multiclass Single Server queueing system, 
we have managed to derive explicit upper bounds on some performance measures such as sojourn times 
and workloads. Then we showed that the bounds implied by the Law of the Iterated Logarithm are appli- 
cable for the underlying stochastic queueing system leading to explicit and non-asymptotic performance 
bounds on the same performance measures. We are not aware of any other method of performance 
analysis which can provide similar performance bounds in queueing model of similar generality. 

We have just scratched the surface of possibilities in this paper and we certainly expect that our 
approach can be strengthened and extended in multiple directions, some of which we outline below. 
First we expect that our approach extends to even more general models, such as, for example multiclass 
queueing networks or more general processing networks [HarOOj . The performance bounds can be 
obtained perhaps again by introducing linear constraints implied by probability laws and using some sort 
of a Lyapunov function for obtaining bounds in the resulting robust optimization type queueing model. 
Another important direction is identifying new probability laws which lead to tighter constraints than the 
ones implied by the LIL. Ideally, one would like to be able to obtain bounds which faithfully represent the 
scaling behavior of the performance measures of interest in the heavy traffic regime as the (bottleneck) 
traffic intensity p converges to the unity. Further, it would be interesting to obtain performance bounds 
on the tail probability of the performance measure of interest, perhaps by constructing constraints 
implied by bounds on the tail probabilities of the underlying stochastic processes. For example, perhaps 
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one can obtain large deviations type bounds by considering the linear constraints implied by the large 
deviations bounds on the underlying stochastic processes. Deeper connection between the results of this 
paper and the results in the adversarial queueing theory and the related queueing literature is worth 
investigating as well. 

Finally, we expect that the philosophy of replacing the probability model with implications of the 
probability model will prove useful in non-queueing contexts as well, whenever one has to deal with the 
issues of stochastic analysis of complicated functionals of stochastic primitives. 
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Appendix. Preliminary technical results 

In this section we establish some preliminary technical results. Using <j) as defined by (|8|), we let 
U(x) = —ax + 2b(j){x) + c for some positive constants a, b, c satisfying 



1995. 



b 

- > e 
a 



2e 



(34) 



Lemma 5. U{x) is strictly concave for x > e e . 



Proof. 




+b( — (lnx) 2 — (xlnlnx) 2 + (lnx) 1 ( )(xlnlnx) a(- |-lnlnx) 






< 



since all three terms on RHS above are negative for x > e e 



□ 
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Lemma 6. Assuming (34\ ) and e e > (c/b) 2 , 

U(x)<0 V x > (186 2 /a 2 )lnln(36/a). 

Proof. Since (18b 2 /a 2 ) In ln(36/a) > e e , throughout the proof we restrict ourselves to the domain x > e e 
Since in addition x > (c/b) 2 , we have b(j)(x) > b^fx > c. In this range —ax + 2bcj)(x) + c < — ax+3b<p(x) = 



— ax + 36\/x In lnx. This quantity is less than zero provided 



x \ \ 36 A 

> — = a. 



V In In x J a 

It is easy to check that xj In lnx is a strictly increasing function with lim x _ s . 00 (x/lnlnx) = oo. Let xq 
be the unique solution of xj In lnx = a 2 on x > e e . We claim that xo < 2a 2 In In a. The assertion of 
the lemma follows from this bound. Let A = 2a 2 In In a. Then 

A _ 2a 2 In In a 

In In .4 ln(2 In a + ln (3) a + In 2) 

2a 2 In In a ( *\ 

> - — — — since In a > hr ' a and In a > In 2 

ln(41no) 

2ct 2 In In a 

> — - — — since ax a > m(b/a) > 2e > 4. 

2m(ma) 

= a 2 . 

This implies xq < A and the proof is complete. □ 
Proposition 2. Under the assumption ^34\j 

sup U(x) < 7(6 2 /a)lnln(6/o) + c. 

x>0 

Proof. Since a > 0, then the supremum in sup 2 . >0 U(x) is achieved. Let x* be any value achieving 
maXa;>o U(x). First suppose < x* < e e . It follows from the definition of <fi in ([5D that (j)(x*) = 1 and 
thus U(x*) = —ax* + 26 + c. Using < x* < e e and assumption (j3lj) . it is straightforward to check that 
U(x*) is indeed upper bounded from above by 7(b 2 /a) lnln(6/a) + c. Next, we consider the case x* = e e , 
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and using the fact that a > 0, we obtain U(x*) < 2b ■ \J e e In ln(e e ) + c. It is again straightforward to 
check that the aforementioned bound is upper bounded from above by 7(b 2 /a) lnln(6/o) + c. 

We now consider the case x* > e e . By Lemma x* is the unique point satisfying ^gjp - =0, if it 
exists. The remainder of the proof is devoted to the final case where we obtain 

dx* Vx* In In x* 

Continuing further, (|35p implies 



\/ x* In In x* 6 / 

a. (36) 



In In x* + j— r a 



Note 



x 



In In x* 



> a 2 



< a since minx* > for x > e e 



2 In In x* In x* 

It is easy to check that x/lnlnx is a strictly increasing function for x > e e and lim x ^ OQ (x / minx) = oo. 
(f34|) implies that there exist unique x m i n and x max satisfying 

^min 2 *^max 9 

— a — — = a 



In In x m i n 2 In In x n 

The monotonicity of x/lnlnx implies x m i n < x* < x max . In order to complete the proof of the 
proposition, we will first state and prove Lemmas [7] and [8l 

Lemma 7. x m i n > a 2 In In a and x max < 4a 2 In In a. 

Proof. Let B\ = a 2 In In a. Then 

B\ a 2 In In a 



lnlni^i In ln(a 2 In In a) 

a 2 In In a 2 

< since in ma > 1 tor a > e 

In ma 

= a 2 . 
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Thus since , f is increasing for x > e e , we have x m i n > -Bi and the first assertion is established. 
Let I?2 = 4a 2 In In a. Then 

B2 4a 2 In In a 

2 In In 5 2 ~~ 21nln(4a 2 In In a) 

4a 2 In In a 



21n(21na + ln (3) a + ln4) 

4a 2 In In a . ^ 

> — — ; since In a > ln w a and In a > In 4 

21n(41na) 

4a 2 In In a 

> — — r- since m a > 2e > 4. 

4 In (In a) 

= a 2 . 



Thus, again since x/lnlnx is increasing for x > e e , then the second assertion follows. 
Lemma [3 and x m i n < x* < x max imply 



a 2 In In a < x* < 4a 2 In In a. 



Lemma 8. \J x max In In x max < 4a In In a. 
Proo/. 



V ^max l n In ^max < \/ (4a 2 In In a) In In (4a 2 In In a) by Lemma [7] 



aV4 In In a a/ In (2 In a + ln^ a + In 4) 



< aV '4 In In aw In (41na) since lna > W 3 ) a and lna > ln(e 2e ) > ln4 



< aV 4 In In aV 2 In In a since lna > 4 
and the lemma follows from the last step. 
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We now complete the proof of Proposition [2j We have 

U{x*) < —ax* + 2bV x* In In x* + c 

^ ax m [ n + 26\/x max In In x max + c since x m ; n ^ x ^ x max 

< — ax m i n + 86a In In a + c by Lemma [8] 

< -oa 2 In In a + 86a In In a + c by Lemma [7] 
= 7(6 2 /a)lnln(6/a) + c. 

□ 
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